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Abstract — We analyze the performance of CDMA signature 
optimization with finite rate feedback. For a particular user, the 
receiver selects a signature vector from a signature codebook to 
avoid the interference from other users, and feeds the correspond- 
ing index back to this user through a finite rate and error-free 
feedback link. We assume the codebook is randomly constructed 
where the entries are independent and isotropically distributed. 
It has been shown that the randomly constructed codebook is 
asymptotically optimal. In this paper, we consider two types of 
signature selection criteria. One is to select the signature vector 
that minimizes the interference from other users. The other one 
is to select the signature vector to match the weakest interference 
directions. By letting the processing gain, number of users and 
feedback bits approach infinity with fixed ratios, we derive the 
exact asymptotic formulas to calculate the average interference 
for both criteria. Our simulations demonstrate the theoretical 
formulas. The analysis can be extended to evaluate the signal- 
to-interference plus noise ratio performance for both match filter 
and linear minimum mean-square error receivers. 

I. Introduction 

In a direct-sequence code-division multiple access (DS- 
CDMA) system, the performance is mainly limited by the 
interference among users. To minimize the interference, every 
particular user wants to select a signature vector from a 
signature codebook to avoid the interference from other users. 
In this paper, we assume that the receiver (base station) has 
the perfect information of the signatures. It selects a signature 
for a particular user according to some criterion, and feeds 
the corresponding index to this user through a feedback link. 
We also assume that the feedback link is error-free and rate 
limited. Due to the finite feedback rate, there is a performance 
degradation compared to the infinite feedback rate case. We are 
interested in quantifying the effect of finite rate feedback. 

This problem has been studied in [1]. A randomly con- 
structed signature codebook is assumed in [1] where the 
codebook entries are independent and isotropically distributed. 
The interference signature matrix is assumed to have indepen- 
dent and identically distributed (i.i.d.) Gaussian elements. A 
particular user chooses the signature vector from the signature 
codebook to maximize signal-to-interference plus noise ratio 
(SINR). For the matched filter receiver, this criterion is equiva- 
lent to select the signature to minimize the interference. In [1], 
an asymptotic lower bound is given on the average interference. 



The main contribution of this paper is to derive the exact 
performance limit. In this paper, we use the average interference 
as a performance measure, which is independent of specific 
receivers and applications. We consider two signature selection 
criteria. One is to minimize the interference from other users, 
same as the one in [1]. The other one is more intuitive. We 
select the signature vector to match the weakest interference 
directions, or equivalently, to be as orthogonal as possible to 
the strong interference directions. To analyze the corresponding 
performance, we let the processing gain, number of users and 
feedback bits approach infinity simultaneously with fixed ratios. 
By asymptotic analysis, we derive lower bounds and upper 
bounds on the average interference for both criteria. For each 
criteria, the asymptotic upper bound meets the asymptotic lower 
bound. Therefore, these bounds provide the exact performance 
limit. The corresponding analysis can be extended to evaluate 
the SINR performance for both match filter and linear minimum 
mean-square error (MMSE) receivers. 

II. System Model 

In a sampled discrete-time symbol-synchronous DS-CDMA 
system, the received vector can be written as 

m 

where Bj e C and Sj E C"^^ are the transmitted symbol and 
the signature vector for user j respectively, and W e i^nxi 
is the additive white Gaussian noise vector with zero mean 
and covariance matrix a^I. The processing gain (length of the 
signature vector) is n, and m is the number of users. We also 
assume that the transmitted symbols Bj's are independent and 
with the same power (variance) 1. 

We assume that the receiver has perfect knowledge about 
the signature vectors s/s. For a particular user, without loss of 
generality, user 1 is assumed, the receiver selects his signature 
to avoid the interference from the other users. It feeds the 
corresponding index back to user 1 through a finite rate and 
error-free feedback link. The rate of the feedback link is 
assumed to be up to Bih bits. In order to accomplish this, a 



signature codebook B with size 2^*^ is declared to both the 
receiver and user 1. 

We assume that the signature codebook B is randomly 
constructed. Specifically, B — {vi , • • • , VjRft }, where v^, = 
Zfe/||zfe||, Zfe = [zi^k,--- ,Zn,fc] and z^^k are i.i.d. C7V(0,1) 
for all 1 < z < n and 1 < fc < 2^^. In this way, 
it is guaranteed that v/c's are independent and isotropically 
distributed unitary complex vectors. It has been shown that the 
randomly constructed codebook is asymptotically optimal [1], 
[2]. 

In this paper, we use the average interference as the per- 
formance measure. Let S e C'^^"-!) be the interference 
matrix for user 1 , whose columns are the interfering signatures 
S27'-- ,Sm- We assume that S has i.i.d. complex Gaussian 
entries with zero mean and variance i, same as the assumption 
in [1]'. For a given interference matrix S, the interference to 
user 1 is defined by 



Theorem 1: Define 
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The average interference is defined by / = Eg [Eg [Is]]- 

In this paper, we consider two types of signature selection 
criteria. The first one is to minimize the interference from other 
users, i.e.. 



Si = arg min v[,SSWfe. 



(1) 



The second one is to select the signature vector to match 
the weakest interference directions (or equivalently, to be as 
orthogonal as possible to the strong interference directions). 
Let d be the multiplicity of the smallest singular value of S. 
Let u„_(j+i, u„_d+2, • ■ ■ I u„ be the d left singular vectors 
of S corresponding to the smallest singular value and \]d = 
[un-d+i ■ ■ ■ u„]. The direction matching criterion is 



si = arg max vj^UdUjjVfc. 



(2) 



For both criteria, we shall derive the asymptotic performance 
limit in Sections Ull] and IIVI respectively. The corresponding 
analysis can be extended to SINR performance evaluation for 
both match filter and linear MMSE receivers [2]. 

III. Analysis for Interference Minimization 

This section is devoted to calculate the average interference 
for the interference minimization criterion in Q. By letting the 
processing gain, number of users and feedback bits approach 
infinity simultaneously with fixed ratios, we derive the exact 
performance limit. The result is given in Theorem ^ 

'it is more natural to assume that the columns in S are independent and 
isotropically distributed unitary complex vectors. However, the asymptotic 
statistics of S are the same for both assumptions. We adopt the assumption in 
[1] for fair comparison. In Section|V] we shall show that the difference between 
these two assumptions is indistinguishable for relatively large systems. 
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for a r > 1, where A* = (1 ± a/t)^- For convenience, define 
X[ ^ \^ \f n < m and A^" = if n > m. For any x £ 
{X^ , A+) and a £ 0, — ^— , define 
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and 



(x, a)= / log(l + a(A-.x))d^A 
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ip (x) = max ij; {x, a) 
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Let 71, TO and i?fb approach infinity simultaneously with fixed 
ratios r = max (n, to) / min (n, to), f = min(n,TO)/ri and 
c = R[b/n. For any < c < oo, there exists an Xc G (A^, A+) 
such that clog 2 — -0 (xc) and 

lim = fxc. (4) 

(n,m,/?fb) — >co 

Remark 1: Theorem is only valid when < c = Rib/n < 
oo. However, it also provides the exact performance limit when 
c — > 0+ or c +00. Elementary computations show that as 
c 0+, fXc ^ f A = f J^- Xd/j-x the average eigenvalue, and 
as c ^ +00, fXc rX~[ the minimum eigenvalue, which are 
consistent with intuition. 

The essential idea behind Theorem [2 is the same as that 
behind the standard large deviation technique. The lower bound 
is derived by Chebyshev's inequality and the upper bound is 
derived by the twisted distribution. Similar to the result in large 
deviation technique, the asymptotic lower and upper bounds are 
identical. We shall outline the proofs for the lower and upper 
bounds in Section UlI-AI and IIII-BI respective Iv. 

As a beginning, we express the average interference in a 
convenient form. For a CDMA system with finite n and m, the 
average interference is given by 



= E.. 
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where J*-"^ is used to emphasize that they are for finite n and 
TO. Let H„ be an n X TO matrix whose entries are i.i.d. complex 
Gaussian CM (0, 1). Obviously, the statistics of S is the same as 
^H„. Therefore, SS^^ = -iH„Hj, where r = min (n, to). 



Let Xi be the i^^ eigenvalue of the matrix -iHnHjj. We have A. The Asymptotic Lower Bound 
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where 



Note 



(a) follows from the random construction of the signature 
codebook B, 

(b) follows from the singular value decomposition of 
iH„Ht„, and 

(c) follows from the fact that z^, and Uzj. are statistically 
equal for any n x n unitary matrix U [3]. 
that given A, the random variables 



Eti^d^^.fclVELil^^l". 1 < < 2«"% are i.i.d.. 
Denote the corresponding conditional distribution function by 
F„ (x |A), then 

It is worthy to keep in mind that this distribution func- 
tion is function of A. Due to the independence of 

Etl^^\^^M^ /Etl\^^M^' 1 < fc < 2«'\ for a givcn A, 
we have 



Pr min 
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Therefore, 
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(1 - F„ (x|A)) 
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where Amin and Amax are the minimum and maximum eigen- 
values of iH„Hjj respectively. 

It is difficult to calculate (|5jl for finite n and to. We let n, 
m and i?fb approach infinity with fixed ratios and derive lower 
and upper bounds on / for large systems. 



The following lemma provides an asymptotic lower bound 
on the average interference. 

Lemma 1: Following the definitions in Theorem ^ let n, m 
and i?fb approach infinity simultaneously with fixed ratios r, f 
and c. For any < c < oo, 

lim > fx^. 

(n,m, /?fb) — *oo 

Due to the length limit, we only sketch the proof. It is based 
on Chebyshev's inequality and the asymptotic behavior of the 
spectrum of a Wishart matrix. Recall that we are dealing with 
a distribution function conditioned on the random vector A, we 
need to define some "good" set of A, say A'^, which will appear 
soon. By Chebyshev's inequality, it can be proved that 

F^{x\\) = Pr (f^(A,;-a;)|z,|' <0|A ) 



< 



diJ,2 



= exp log (1 + g {Xj - x))j 

for Va e (o, ^ j and A e A'^, where the set is defined 
by 

Al ^ {A: |^„(A,x,a)-V(a^,a)| <ei} 

n{A: I Amin 1 <e2}, (6) 

1 " 

ipn (A, x, a) = - log (1 + a {Xi - x)) , 
n ^-^ 

and the positive numbers ei and 62 are small enough. By the 
asymptotic behavior of the spectrum of a Wishart matrix, it can 
be shown that Pr (A^J;) 1. 

Now take a small e > such that Xc— e > A^. Let x ~ x^—t. 
It can be proved that we can always find an a e ^0, — 
such that i\j {x, a) > clog 2. Then for V(5 > 0, 



(l-F„(x|A)) 
>(l- 

p (^2™ log (l - e""" ^ log(l+a(A,-i:))^ 
> exp (^_e-«['A(:';,a)-^l-clog2] ^^^Q (^^^^ 
(b) 



(a) 



(7) 



on for n large enough, where 

(a) follows by Taylor series expansion, and 

(b) follows from the fact that we are able to choose ei > 
small enough such that tp {x, a) — ei — clog 2 > 0. 



Then 



> Ea 



> Ea 



{1-Fr,{x\\)f dx 



(1 - K {x\\)f dx, Al 



(1 - 5) dx, Al 

= {l-5) {x,-\i -e-e2)-V,[Al) 
-f {1-5) (xc - - e - £2) , 



where 
(c) 



follows by reducing the integration domain of a non- 
negative function, and 
(d) follows from Q and the fact that F„ {x\\) is a non- 
decreasing function in x. 
Therefore, 



lim/(") 



r lim Ea 



{l-Fn{x\X)Y dx 



By taking 6, e and e2 arbitrarily small, we have lim/'") > fxc- 
B. The Asymptotic Upper Bound 

For the interference minimization criterion in Q, the asymp- 
totic upper bound on the average interference is given in 
Lemma |2] 

Lemma 2: Following the definitions in Theorem ^ let n, m 
and i?fb approach infinity simultaneously with fixed ratios t, f 
and c. For any < c < oo, 

lim < fXc. 

(n,m,i?.fb) — *oo 

Due to the length limit, we omit the detailed proof. A sketch 
of the proof is given in the below. 

To prove the upper bound, roughly speaking, it is sufficient 
to show that for Ve > and V(5 > 0, if n is large enough, we 

2CTI 

can upper bound (1 — F„ (a::|A)) uniformly by 5, i.e., 

(1 - Fn (a;|A))^ < 5 for all x> Xc + e, 
on the "good" set A\ (|6j. Since Pr {A^) 1, 



(8) 
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(1 - F„ {x\X) f dx 



< rlEx [Xnun] + Ea 



Idx 



+Ea / ldx,nx-Al 

+Ea / " 5dx,Al 

f {Ea [x, + e] + (A„,ax - - e) (1 - Pr (A^)) 
+S {Xn... - X, - e) Py{AI)} 

f{Xc + € + S (Ainax - Xc ~ e)} ■ 



By taking e and S arbitrarily small, we have 

lim/(") < fXc- 

The essential tool used to prove ^ is the twisted distribu- 
tion [4]. This tool is the main tool in proving the lower bound of 
Cramer's theorem [4], a basic result of large deviations. How- 
ever, there is a fundamental difference between the standard 
large deviation technique and our approach. While in Cramer's 
theorem one considers the sums of i.i.d. random variables, 
here we consider the sum of (Ai — x) \zi\^, where the random 
variables are conditional independent but not identically dis- 
tributed and the condition itself is a random vector. While the 
conditional independence requires us to discuss the statistics of 
J2 {Xi — x) \zi\^ on the set A^, the non-identical distribution 
brings the major difficulty. That is, the twisted distribution may 
or may not be well-defined. To overcome this difficulty, we 
have to discuss two types of x and define two types of twisted 
distributions respectively. 

We define two types of x and two types of twisted distribu- 
tions as follows. Let be the a such that (x) — ip {x, ax)- 
The set of x of the first type is defined by 



Xi = jxe (Ar,A+) 



e 0, 



1 



X- Xt 

The set of x of the second type is defined by 

1 



X2 



e (At , A+) ; ax 



a; 



It can be proved that the Xc in Theorem Q] is either in Xi or in 



A2. If X e Xi, then ax < 



and E, 



well defined on the set A^ with £2 small enough. Then we are 
able to define a twisted distribution measure 



dflz 



-aY:i\-^)\z,\^ 



Ez [e 



where d/i^ is the probability measure for the random vector z. 
However, if x € X2, E^ e-'^T.i^i-x)\^i\^ is not well defined 
on the set no matter now small £2 > we choose. For this 
case, we have to define the twisted distribution in a "truncated" 
way. Define the Af -truncated measure for z as 



df^^' = '[[^z,elo,Ai]dpz 



Then the A/-truncated twisted distribution measure is defined 
by 

g-ai:(A.-a;)|z.|^ 



r.M A 



E„ 



-dp^ 



M 



It can be verified that dp^^ is always well defined on the set 
Al for a finite M > 0. In the proof of (|8}, we need to choose 
an M sufficiently large. 

With the twisted distributions, (|8} can be proved. Here, we 
only outline the proof for Xc E Xi, the simpler case. Assume 
that Xc G Xi. For an e > small enough, let a; = + ■§. It 



can be proved that x E Xi and i/j (x) = ij;{x,ax) < clog 2. 



For V(5i > and a y > x 



it can be proved that 



Fn (y|A) > Pr (J2 (A, - x) \z,\^ < nea) (1 - 6^) 

on the set with small enough 63 and large enough n. But 
for any given €3 > 0, a further lower bound can be derived as 
follows. 

Pr (J2 {>^^-x)\z,f <ne3^ 

{Xi ~ x) \z^f e n{-62,S2)j 



(a) 

> 

(b) 
> 



Pr 



-n52 



> 



exp {-n [S' + ijj (x, a^)]} P^^ {B^ 



on the set A\ with large enough n, where 

(a) holds by choosing 82 < 

(b) holds by defining 

= |z : -a^ ^ {h - x) \zi\^ e n {-62, 62)^ , 

(c) follows from the definition of the twisted distribution 
for X E Xi, and 

(d) follows by letting 6' = t\ + i52, where (Pz) is the 
probability of under the twisted distribution. 

We want to calculate P^^ (-Bz)- By studying the asymptotic 



behavior of log 



it can be proved 



that for V(5" > 0, P^,^ (B^) >\-8" on the set with large 
enough n. Therefore, for V(5"' > 0, 

Fn (y|A) > e-"['^(^)+'''] (1 - 6"') 

on the set with large enough n. Now we choose b' small 
enough such that ^/i (a;) + 28' < clog 2. Then it can be proved 
that, for > 0, 



< 
< 



[l-Fn{y\X)f 

exp (1 - <5"') e-"['^(^)+*'-=i°g2] 
6 



on the set A'^ with large enough n. Without loss of generality. 



we take y = Xc 



e > X — Xr 



Since F„ (2/|A) is a non- 



decreasing function, we have uniform boundedness, 
[1 - F„ (y|A)]^'" < S for al\y>Xc + e 
on the set with large enough n. This is (|8}, what we want. 

IV. Direction Matching Criterion 

In this section, we shall analyze the performance correspond- 
ing to the direction matching criterion in Again, by letting 
n, m and i?fb approach infinity simultaneously with fixed ratios, 
we derive the exact performance limit. The result is given in 
Theorem |2l 



Theorem 2: Following the definitions in Theorem ^ let n, 
m and i?fb approach infinity simultaneously with fixed ratios 
r, f and c. For any < c < 00, 



lim = 

(n,m,i?fb) — >OQ 



Ar(i 



2-^) + A2- 



if n < m 
if n > TO 



(9) 

where X — Xc < ^ satisfies D {/ip \\ fJ'xJ = clog2 and 
D{^J'f II Majc) - ^logi; + (l-'^)logTf^ is known as the 
relative entropy. 

Remark 2: Elementary calculations show that the asymptotic 
average interference lim /, as a function of c, converges to 
the average eigenvalue and the minimum eigenvalue as c — > 
0+ and c ^ 00 respectively. These results are consistent with 
intuitions. 

The proof of Theorem |2l is based on the observation that 

n—d 

/(") = f ^ Eu [Eb piu,ujsi] ] Ea [A,] 



+f 
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i=^n—d-\-l 
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E, 



s|u,jU^Si 



Ea [A„ 



where Ai > • • • > A„ are the singular values of iH„Hjj and 
Uj is the singular vector corresponding A^. For n < m (full 
rank) case, d = 1 with probability 1. We select the signature 



s|u„uJjSi 



can 



Si to match u„. The corresponding Eu 
be calculated based on our previous resiilts in the Grassfnann 
manifold [5]. For n > m (deficient rank) case, d — n — m 
with probability 1. We need to choose the signature Si to 
match the plane generated by = [un-d+i ■ ■ ■ u„]. By large 



deviation technique, the corresponding Eu 



Eg 



can be evaluated. The detailed derivation is given in [2]. 

V. Simulations and Discussion 

Fig. n shows the simulation results to demonstrate the 
asymptotic performance formulas and (|9jl for both criteria. 
Fig- E^) ^b) are for n < m (full rank) and n > m 
(deficient rank) cases respectively. From the simulations, we 
can observe that the simulated average interference for both 
criteria (x markers for interference minimization and circles 
for direction matching) converges to the asymptotic results (the 
solid line for interference minimization and the dashed line for 
direction matching) as n, ni and i?fb approach infinity with 
fixed ratios. Simulations also show that direction matching is a 
sub-optimal criterion. 

We also compare our formula with the bound in [1], denoted 
as SH bound in Fig^ In [1], an asymptotic lower bound on the 
average interference is given for the interference minimization 
criterion. It is plotted as the dotted line in Fig.^ Note that when 
T = 2 and n < to, with infinite feedback rate (c = 00), 
should converge to the minimum eigenvalue (l — \/2) w 0.17. 
The bound in [1] is below this value even when c is relatively 
small (c > 2.5 in FigQ. Generally speaking, the bound in [1] 
under-estimates the interference while our asymptotic formula 
(Hi gives the exact performance limit. 
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(b) Deficient Rank (n > m) Case 
Fig. 1. Simulations for both signature selection criteria 
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Fig. 2. Comparison of two types of interference matrices 



As mentioned in the system model section, we assume that 
the interference matrix S has i.i.d. complex Gaussian entries 
with zero mean and variance i for fair comparison, while 
it is more natural to assume that S has independent and 
isotropically distributed unitary complex columns. Fig |2] gives 
the difference between these two statistical assumptions, where 
the simulations are based on the interference minimization 
criterion. For small n and m, these two different assumptions 
give two different results. However, as n and m increase, 
for example, n = 16 and m = 32, the difference becomes 
indistinguishable. Indeed, the asymptotic statistics of these two 
types of random matrices are identical. The results Q and 
are the exact asymptotic performance for both interference 
statistical assumptions. 

VI. Conclusion 

In this paper, we quantify the average interference as a 
function of finite feedback rate for CDMA signature optimiza- 
tion problem. Two signature selection criteria, i.e., interfer- 
ence minimization and direction matching, are analyzed. By 
letting the processing gain, number of users and feedback 
bits approach infinity with fixed ratios, we derive the exact 
asymptotic formulas to calculate the average interference for 
both criteria respectively. The asymptotic results are valid for 
both the Gaussian interference matrix and the interference 
matrix with independent and isotropically distributed columns. 
Furthermore, the corresponding analysis can be extended to 
SINR performance evaluation for both match filter and hnear 
MMSE receivers. 
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